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We report on the Lattice determination of the semileptonic form factors by lattice QCD. Comparison with the light-cone QCD sum rules 
are made for B — * nlv, B — > p/v semileptonic decays. 



1 Introduction 

Computation of the semileptonic form factors for B — > 
n(p)lv, B -> D w /v and B -> K^l + t decays is a key to 
the determination of the CKM matrix elements |V u b|, |V c bl 
and | Vts I ■ In this report, I review the recent results by lat- 
tice QCD methods. I will mainly focus on the semileptonic 
decay B — > nlv, because lattice calculations in various ap- 
proaches are already at hand [ Q] |2 El EJ and because 
most of the problems which exists in this decay are com- 
mon to other processes. 

Lattice QCD calculations of B — > nlv form factors suffer 
from (1) the error from the large heavy quark mass mq in 
lattice unit, (2) the chiral extrapolations in the light quark 
mass m q using the simulation results with relatively large 
mass region m q > m s /2, and (3) the limitations from the 
accessible kinematic range of q 2 from statistical and sys- 
tematic errors. 

In order to solve the first problem, two different approaches 
are made. One is to avoid the large discretization errors of 
O(amQ) by computing the form factors with the conven- 
tional relativistic quark action for charm quark mass re- 
gion, and then extrapolate the results in 1 /rriQ. The other is 
to use HQET effective theory with 1/niQ corrections. Since 
both approaches have their own advantages and disadvan- 
tages, it would be important to have both results and study 
whether they give consistent results within quoted errors. I 
will present some of the major calculations from these two 
approaches and discuss the consistency of the results. 

The second problem already gives a source of errors in the 
quenched calculations but it would become even more se- 
rious in unquenched QCD. Unfortunately, there are still no 
results in the unquenched QCD. Since the chiral limit of 
the form factors is ill-defined in the quenched QCD, all we 
can do with the present lattice data is to discuss the light 
quark mass dependence in the intermediate mass regime in 
the quenched QCD. However, one may be able to give an 
estimate of the low energy coefficients of the chiral per- 
turbation theory based on the present lattice data. I will 
briefly review new studies of quenched chiral perturbation 
theory (QChPT) and partially quenched chiral perturbation 
theory (PQChPT) [ |5) for Z? — > nlv form factors, which 



give a phenomenological estimate on the quenching errors 
and chiral extrapolations. These studies will be even more 
useful once new calculations in unquenched QCD will be 
made. 

The third problem arises because in semileptonic decays 
large energies are released to the final states so that the lep- 
ton pair invariant mass q 2 can range from to (ms - m^) 2 . 
However, due to the discretization errors of O(aE) as well 
as the statistical errors which grow as ~ exp(const x (E - 
m^)i) where E is the energy of the pion, lattice QCD can 
cover only large q 2 region. The dispersion relation is a pos- 
sible solution to give bounds for smaller q 2 region [ 7 1 . The 
light-cone QCD sum rule (LCSR) predicts form factors for 
small q 2 , which is complementary to the lattice results. I 
will give a comparison of the recent LCSR results with the 
lattice results to see whether they will give consistent re- 
sults. 

1 also review the form factors in other processes. Some 
of the recent work on the lattice QCD calculations of the 
B — > plv form factors in relativistic formalism are pre- 
sented. Very precise calculations of semileptonic form fac- 
tors for B — > D w /v at zero recoil and the calculations of 
the slope of the Isgur Wise function are presented. 

2 B^nlv 

The exclusive semileptonic decay B — > nlv determines the 
CKM matrix element \ V u b\ through the following formula, 



If = ^i 2 i^ivv)i 2 , 

where the form factor / + is defined as 



(1) 



(n(k n )\qy»b\B(p B )) 



/V) 



(Pb + k K f 



2 2 

mt - mi 



/V) 



q 2 



(2) 



with q = p B - k n and q 2 = m 2 B + m 2 - 2m B v ■ k„. The 
following parameterization proposed by Burdman et al. [ 
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(n(k n )\qfb\B( V )) = 2 



/i(v • k n y + f 2 (v K)—, 

v ■ k n 



(3) 



is also useful for discussing the heavy quark symmetry and 
the chiral symmetry of the form factors in a transparent 
way. 

2.1 Lattice results 

Lattice calculation is possible only in limited situations. 
Spatial momenta must be much smaller than the cutoff, i.e. 
\Pb\, \kn\ < 1 GeV. This means v-k„ = E„ < 1 GeV orequiv- 
alently q 2 > 18 GeV 2 . Another limitation is that due to 
the slowing down, simulations with very small light quark 
masses are difficult so that usual mass range for the light 
quark masses in practical simulations is m s /3 < m q < m s 
or m„ = 0.4 ~ 0.8 GeV. Therefore in order to obtain phys- 
ical results chiral extrapolations in the light quark masses 
are necessary. 

So far all the lattice calculations of the form factors are 
done only in quenched approximation. APE collaboration 
[ [2 and UKQCD collaboration computed B — > nlv form 
factors for a fine lattice with the inverse lattice spacing 
aT l ~ 2.7 GeV. They used relativistic formalism for the 
heavy quark and extrapolated the results of heavy-light me- 
son around charm quark masses to the bottom quark mass. 
Fermilab collaboration [ 3| used the Fermilab formalism 
for the heavy quark and computed the form factors on three 
lattices with or 1 = 1.2 ~ 2.6 GeV . JLQCD collaboration 
[|4) computed the form factors using NRQCD formalism 
for the heavy quark on a a~ l = 1.64 GeV. NRQCD col- 
laboration [ 5 1 also used NRQCD formalism for the heavy 
quark and an improved light quark action (D234 action) 
on a anisotropic lattice with aT x = 1.2 GeV (spatial), 3.3 
GeV (temporal). In all of these calculations the light pseu- 
doscalar meson masses are 0.4 ~ 0.8 GeV. Fig. 1 shows 
the result by different lattice groups, f + {q 2 ) agrees within 
systematic errors while f l> (q 2 ) shows deviations among dif- 
ferent methods. 

The reason for the discrepancies in f° can be attributed to 
the systematic error in the chiral extrapolation and heavy 
quark mass extrapolation (interpolation) error. In the fol- 
lowing, we examine these errors in more detail. Light 
quark mass m q dependence of form factors with fixed spa- 
tial momenta ap = ||(1,0, 0) is shown in Fig. 2. In con- 
trast to the JLQCD data, Fermilab data shows a significant 
increase towards the chiral limit. Large difference in Fer- 
milab results and JLQCD results for f° in the chiral limit 
arises from different m q dependence, but the raw data for 
similar quark masses are not so different. Shigemitsu et al. 
studied the mass dependence of f\ + fz and find similar be- 
havior as JLQCD. Further studies to clarify the light quark 
mass dependence are required. 



a; 

+ 



o.o 



UKQCD 




APE 


i 


□ Fermilab 




JLQCD J 




- NRQCD i *JK 


/ 


LCSR ifJfff 




- LCSR ^jffl d 











10 , 15 , 20 25 30 

q (GeV ) 



Figure 1. B — > nlv form factors by different lattice groups. 

Fig. 3 shows I 1Mb dependence of form factors <D> + = 
Vmi/ + at v • k„ = 0.845 GeV for APE and JLQCD col- 
laboration data. It is found that the difference of APE 
(UKQCD) vs JLQCD (NRQCD) for/ arises from the ex- 
trapolation in 1/M. Linear extrapolation in 1/M is consis- 
tent, while the quadratic extrapolation gives higher value. 
The quadratic extrapolation 1/M is chosen for APE's result, 
since higher value gives better agreement with the soft pion 
theorem. Simulations with static heavy quark may resolve 
the problem. 

The error of the form factors in the present calculations is 
around 20%. Some of the major errors are the quenching 
error, chiral extrapolation error statistical error in all cal- 
culations. In addition, a large discretization error appears 
in JLQCD results and a large 1/M extrapolation error is 
contained in APE and UKQCD results. 

There are several proposals to improve the form factor de- 
termination. The quenching error can be resolved only 
by performing the unquenched calculations. Recently, 
JLQCD and UKQCD collaborations has accumulated «y ■ - 
2 unquenched lattice configurations with <9(a)-improved 
Wilson fermions and nj = 2 + 1 unquenched configurations 
with improved staggered fermions have been produced by 
the MILC collaboration. These unquenched QCD data 
should be applied to form factor calculations. 

In order to reduce the chiral extrapolation error, simulation 
with even smaller light quark masses are necessary. For 
Wilson type fermions, simulations with m„ < 0.4 GeV will 
be very slow and also appearance of exceptional configura- 
tion may prevent the simulation for very light quark mass 
range. On the other hand, MILC collaboration is now car- 
rying out simulations with m„ - 0.3 - 0.5 GeV, which cor- 
responds to m q = l/5m s -l/2m s Since rtf = 2+1 sim- 
ulations are performed by taking the square root or quar- 
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Figure 2. The left triangle symbols denote the Fermilab results, 
while the circle, square, and diamond symbols denote the JLQCD 
results for different heavy quark masses. A clear rise towards the 
chiral limit is observed for Fermilab data. 



tic root of the Dirac operator, one possible concern might 
be the flavor symmetry breaking effect which still exists 
in the improved stagger fermions. It would be important 
to understand how this flavor breaking affects distorts the 
chiral behavior. Although it is possible to simulate with 
very light quark mass without encountering any theoretical 
problems, the simulation cost at this stage is very high. A 
breakthrough in the algorithm is needed. 

Using the heavy quark symmetry is another way for im- 
provements. Since CLEO-c experiment can measure form 
factors for D — > nlv to a few percent accuracy, their results 
will be a good approximation for the B — > nlv form factors. 
Then the task for lattice QCD is to predict the 1/toq de- 

dT(B -> nlv) 

pendence of the form factors. The B/D ratio 

dT(D — > nlv) 

with the same recoil energy v ■ k„ would be a nice quantity 
to measure on the lattice, since a large part of the statisti- 
cal error, the perturbative error and the chiral extrapolation 
errors are expected to cancel in this ratio . 

2.2 QChPT and (PQ)ChPT 

Becirevic et al. [ |6) made estimates of the quenching ef- 
fect and the chiral extrapolation error from low energy ef- 
fective meson theory, namely QChPT and (PQ)ChPT for 
quenched QCD and full QCD respectively. 

In their analysis the non perturbative low energy coupling 
constant are (1) a, g and / for fs B*Bn coupling and f„, (2) 
L4 and L5 for the two 0(p 4 ) terms, and (3) mo and g' for the 
two parameters of quenching effect. The parameters are es- 
timated by collecting all the knowledge from the quenched 
lattice QCD simulation, the full QCD simulation, the ex- 
perimental values and large N order estimation. 



Figure 3. 1/M scaling with fixed v ■ k n = 0.845 GeV. The circle 
symbol denotes the JLQCD data, while the diamond symbols de- 
notes the APE data. Upward and downward triangles correspond 
the extrapolated results of APE data by linear and quadratic func- 
tions in 1/M. 

They found that the quenching errors Q2 = (f( ul1 - 

^uench )/f full ^ Q ^ _ + fry/nil ^ +/ 2 )*^))/(/ 1 + 

fi) are as large as 25%-50% depending on v ■ k n and g' . 

They also estimated the chiral extrapolation error by com- 
paring the result with linear extrapolation and linear plus 
chiral log corrections for unquenched theory. They found 
that the linear extrapolation of m n > 0.4 GeV data gives an 
overestimate of f ± oc f 2 by 2-7% for v ■ k n — 0.19 GeV and 
5-15% for v ■ K = 0.54 GeV. 

They proposed an extrapolation strategy in the partially 
quenched theory for m s n ea > 0.4 GeV simulation data in 
which one first takes a linear extrapolation of m sea with 
fixed nivaience and then make an extrapolation in m va i ence 
with a linear + log form. In this method they estimate that 
the chiral extrapolation error in unquenched theory in un- 
der 5% level. 

2.3 q 2 dependence from the dispersive bound and 
LCSR 

Model independent bounds for the whole q 2 range can be 
obtained with the dispersion relation, perturbative QCD, 
and the lattice QCD data [[7). Fig. 4 shows the pioneer- 
ing result by Lellouch [ 7 |. Reducing the lattice errors or 
having other inputs would significantly improve the results. 
More elaborate studies along this line would be important. 

The light-cone QCD sum rule (LCSR) \ MD1 [TT1 ITU can 
also give form factors for small q 2 region. The basic idea 
is to compute the following matrix element CF V 



CF V = 



"I 



d 4 ye l<1 >\n(p)\T[qy,b](y)j B (0)\0}, 



(4) 



4 
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the correction from higher twists ( twist 3) is ~ 30%. The 
results for q 2 < 14 GeV 2 are well fitted by 



/V) = 



F(0) 



1 — aq 2 /m 2 + b(q 2 /m B ) 



(7) 




Light-cone QCD sum rule results for q 2 < 14 GeV 2 can 
also be fitted by the pole dominance ansatz. 



/V) 



1 - q 2 /m 2 B , 



(8) 



where c = fB-gBB'n/(2mB-) - 0.414^gg[g plus systematic 
errors. 



Fig- EJ] shows the result by the light-cone QCD sum rule. 
It is remarkable that the light-cone QCD sum rule give con- 
sistent results with lattice QCD. 



3 B^plv 



Figure 4. Dispersive bound for f° and / + with UKQCD lattice 
data. Model-independent QCD bounds with 90%, 70%, 50% and 
30% confidence levels are given by the pair of curves. Figure 
taken from Ref. r|7|. 



in two different methods and equate the results. The two 
computational methods are : (1) light cone expansion 
which is expressed by the pion light-cone wavefunction 
4> K {u) and (2) the dispersion relation which takes the fol- 
lowing sum of the physical poles 



CF V 



>n 2 B fB 



/V)" 



1 



pI 



+ higher poles, 



(5) 



where higher poles are suppressed by Borel transformation 
with M and approximated by the light-cone expansion re- 
sults above a threshold si. 

The theoretical input parameters are the parameters a,'s 
of light-cone wavefunctions in the Gegenbauer polynomial 
expansion, 



6m(1 



m)[1 + ci2CV 2 (2u ■ 



!) + ••■], 



(6) 



the B meson decay constant /g, the b quark mass nib, the 
threshold s 2 r and the parameter M for the Borel transfor- 
mation. 

We here give the new results of Ref. Ml 21 as an example. 
It is found that the radiative correction is about 10% and 



Recently, UKQCD collaboration [ 1 3 ] and SPQcdR col- 
laboration [ 14 1 started studies of B — > plv form factors. 
Both collaborations use (9(a)-improved Wilson action for 
the heavy quark and extrapolate the numerical results of 
mg ~ m c towards the physical b quark mass. The lattice 
spacings are a = 2.0 and 2.7 GeV for UKQCD and a 
2.7 and 3.7 GeV for SPQcdR. 

UKQCD fits the lattice data for q 2 > 14 GeV 2 to the fol- 
lowing form 



i dr 

W^dq 1 



G 2 F q 2 [A(q 2 )] 1 / 2 

— 1 — (" + bir - O). 



1927T 3 /. 



The fit coefficients are a = 38+° ± 5 GeV and b = + 2 + 
1, where the first error is statistical and the second is the 
extrapolation error for both a and b . 



SPQcdR collaboration obtains form factors for q > 10 
GeV 2 . They find the results which is consistent with the 
light-cone QCD sum rule results. 



4 B D { *Hv 

One can extract \V c /,\ from the B 
decay near zero recoil as 



dr . 

da) 



D^lv semileptonic 



(9) 



where co = v • v' and f^ D (.i (f are the linear combinations 
of form factors h ± ,h An} . One important outcome from the 
heavy quark symmetry is that the form factor Tb^d^Hv is 
equal to unity at zero recoil up to perturbatively calculable 
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Figure 5. A comparison of B —* plv form factors from lattices 
for q 2 > 10 GeV 2 with the light-cone sum rule (LCSR) results 
q 2 < 10 GeV 2 region. Figure taken from Ref. f l 141 . 



factors and l/m;,, c corrections only starts from the second 
order as, 



1 - 



21a 



If 



(2m c ) 2 2m c 2nih (2nib) 2 



2m c 2nib 



(10) 



Hashimoto et al. N15I computed the double ratio in lattice 
QCD to obtain the coefficients IvJaAp 



{D\cy°b\B){B\by°c\D) 



= I 



1 



1 



{D\cy°c\D){B\by Q b\B) 

(D'\cy°b\B')(B'\by a c\D') _ 
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1 - l P { ) 

2m c 2mi, 



1 - h{^~ - 



1 



2m c 2mi, 



1 



1 



2m c 2mi, 



(11) 



(12) 



I 2 (13) 



The deviation from unity in the double ratio only start from 
1/m 2 including systematic errors. The B — > D* result is 



n. Qi O+0.024 , rv 01 fi+0.003+0.000+0.006 
u - ylJ -0.017 ± u - uiD -0.014-0.016-0.014' 



where the first error is statistical and the second, the third, 
the fourth and the fifth are the errors from the perturbation, 
the discretization, the chiral extrapolation, and quenching. 

The B — > D result is 



Tb^d(D = 1.058 + 0.016 ± 0.003 



where the first is statistical and the second is the error from 
m.b, c and the third is the perturbative error. 

UKQCD collaboration [ 16 1 evaluated the co dependence 
of the form factors B — > D^lv from the 3 point functions 
in quenched QCD. Parameterizing the form factor as 



hj(aj) 



m = i 



-p 2 (oj ■ 



1) + 0((cj 



I) 2 ) 



They obtain the slope 



^ 2 - n on+0.15+0.24 



(14) 



(15) 



where the first and the second errors are statistical and sys- 
tematic, respectively. 

UKQCD collaboration [ 17 1 also computed the form fac- 
tors Kb — > A c lv . Gottlieb and Tamhankar H18I also made 
an exploratory study on he same form factors and obtained 
t;{u>) for finite heavy quark mass. Both groups found that 
the heavy quark mass dependence is small. 

5 Summary 

B — > tiIv form factors are computed with 20% error for 
q 2 > 18 GeV 2 . Different approaches show good agree- 
ments in the form factors for the same mass range of the 
heavy and the light quarks, although the extrapolation in m q 
and 1 /Mq gives deviations in f°. Therefore, a better under- 
standing of the quark mass dependence is required. Studies 
in QChPT and PQChPT suggests that the unquenching er- 
ror may be quite large so that lattice results in unquenched 
QCD are required. Once unquenched lattice results will be 
available, studies in PQChPT suggest that chiral extrapo- 
lation errors can be controlled below 5% level. The form 
factor for smaller q 2 range can be obtained either by the 
dispersive bounds or Light-Cone QCD sum rule with the 
help of nonperturbative inputs. B — > plv form factors are 
in progress. It seems that the form factors in the range q 2 > 
14 GeV 2 are feasible. B — > D^*Hv results are established in 
quenched QCD, whose extension to the unquenched calcu- 
lation should be straightforward. 
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